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AMENDMENT 

Amendments to the Claims 

The claims are amended as shown on the following pages under the heading 
LIST OF CURRENT CLAIMS. The list shows the status of all claims presently in the 
application including any current amendments. This list of claims is intended to 
supersede all prior versions of the claims in the application. Any cancellation of 
claims is made without prejudice or disclaimer. 
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LIST OF CURRENT CLAIMS 

1 . (Currently Amended) A statistical facial feature extraction method, 
comprising: 

a first procedure for creating a statistical face shape model based on a plurality 
of training face images, including: 

an image selecting step, to select N training face images; 

a feature labeling step, to respectively label feature points located in n 
different blocks of the training face images to define corresponding shape vectors of 
the training face images; 

an aligning step, to align each shape vector with a reference shape vector to 
thus obtain aligned shape vectors; and 

a statistical face shape model computing step, to use a principal component 
analysis (PCA) process to compute a plurality of principal components based on the 
aligned shape vectors to create a statistical face shape model, wherein the statistical 
face shape model represents the shape vectors by combining a plurality of projection 
coefficients, and the statistical face shape model computing step includes: computing 
a mean value of the feature points of the aligned shape vectors to define a mean shape 
- 1 » 

_ x = — 2X 

vector x as j_j a 1 where x a is an aligned shape vector , subtracting each 

aligned shape vector x a by the mean shape vector x to form a matrix 
A = [d X) ,d X2 ,...,d XN j and d x = x a — X , computing a covariance matrix C of the 
matrix A, and computing a plurality of principal components according to 
eigenvectors v,. , which are derived from Cv* — ^k V k w * m eigenvalues X\ 

corresponding to the covariance matrix C formed as C = AA T for 1 =k=m where m 
is the dimension of the covariance matrix C for A, 5 , > X, s 2 > ... > ?t s mi to thereby form 

the statistical face shape model; and 

a second procedure for extracting a plurality of facial features from a test face 
image, including: 
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an initial guessing step, to guess initial positions of n test feature points located 
in the test face image, wherein the initial position of each test feature point is a mean 
value of the feature points of the aligned shape vectors; 

a search range defining step, to define n search ranges in the test face image, 
based on the initial position of each test feature point, wherein each search range 
corresponds to a different block; 

a candidate feature point labeling step, to label a plurality of candidate feature 
points for each search range; 

a test shape vector forming step, to do combination of the candidate feature 
points in different search ranges in order to form a plurality of test shape vectors; and 

a determining step, to match the test shape vectors respectively to both the 
mean value and the principal components for computing a similarity, and to 
accordingly assign one feature point corresponding to one, having the best similarity, 
of the test shape vectors as facial features of the test face image. 

2. (Original) The method as claimed in claim 1, wherein in the feature labeling 
step of the first procedure, the feature points are coordinates for corners of eyes and 
mouth on each training face image. 

3. (Original) The method as claimed in claim 1, wherein the feature labeling 
step of the first procedure manually labels the feature points of each training face 
image. 

4. (Original) The method as claimed in claim 1 , wherein the reference shape 
ector is one of the shape vectors. 

5. (Original) The method as claimed in claim 1, wherein the aligning step of 
the first procedure uses a 2D scaled rigid transform algorithm to align each shape 
vector with the reference shape vector. 
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6. (Original) The method as claimed in claim 5, wherein the aligning step of 
the first procedure further comprises the steps of: 

selecting the reference shape vector as x ; = (x il ,y jl ,...,x in ,y in ) and one of 

the shape vectors as x. =(x. 1 ,y jl ,...,x jD ,y jn ); 

computing a squared Euclidean distance E between the vectors Xj and xj based 
on the following equation 

E = (x. - M <N) (a,e)[ Xj ] - t) T (x f - M (N) (a,e)[xJ- t), where 
M (N) (a,0)[x j ] — t is a geometric transform function defining with a plurality of 
transfer parameters to align the vector x j5 the transfer parameters include a rotating 
angle 0 , a scaling factor a , and a shifting vector represented by t = (t x , t y ) , and 

M (N) (a,0) is a 2nx2n diagonal blocked matrix as well as 



M(a,8) 

_> 

M(a,9) = 



acos9x jk -asin9y jk 
asin0x jk +acos0y jk 



for 1 ^k^n, as 



[ acos0 - ctsinG^ 
i^asinB acosGj' 
finding the smallest squared Euclidean distance and corresponding rotating 
angle Q j , scaling factor CL i and shifting vector represented by tj = (t xj , t yj ) to align 

the shape vector xj so similar as the reference shape vector X;; 

computing a sum of smallest squared Euclidean distances after the N shape 
vectors are all aligned so similar as the reference shape vector, ending the aligning 
step when the sum is smaller than a predetermined threshold; 

computing a mean value of the feature points in each block for the aligned 
shape vectors to define a mean shape vector for each aligned shape vector as 
1 N 

x = — X x » wherein x a is the aligned shape vector; and 

N a=l 

assigning the mean shape vector as the reference shape vector and the aligned 
shape vectors as the shape vectors and then repeating the aligning step until all shape 
vectors are aligned. 
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7. (Original) The method as claimed in claim 6, wherein the transfer 
parameters is obtained by a least square algorithm. 

8. (Original) The method as claimed in claim 1, wherein the statistical face 
shape model is a point distribution model (PDM). 

9. (Canceled) 

10. (Original) The method as claimed in claim 1, wherein each shape vector Xj 
consists of n feature vectors Sj located in different blocks, so an average value as 

t = — of the feature vectors s, corresponding to special blocks of all shape 

vectors xj is defined as a feature template. 

11. (Original) The method as claimed in claim 1, wherein in the initial 
guessing step of the second procedure, scaling of initial guess shapes formed by the 
test feature points is aligned so similar as the test face image. 

12. (Original) The method as claimed in claim 10, wherein the candidate 
feature point labeling step of the second procedure further comprises the steps of: 

k 

labeling a plurality of reference points derived from I. = t + £ b p . 

i=. J J 

respectively in each search range, where t is the feature template of block 
corresponding to a search range, pj is j-th principal component of the statistical face 
shape model computed from the training feature vectors, and bj is an associated 
projection coefficient; 

k 

using e =|| I. — t — Zbjpj || 2 to compute an error value between one of the 

j=i 

reference points and the corresponding principal component pj and projection 
coefficient bj; and 
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selecting preceding k smallest error values and defining the k smallest error 
values as feature points of the search range. 

13. (Original) The method as claimed in claim 12, wherein the test shape 
vector forming step of the second procedure does combination of the candidate 
feature points in different search ranges to thus form k" test shape vectors. 

14. (Original) The method as claimed in claim 10, wherein the determining 
step of the second procedure further comprises the steps of: 

using the average value of the aligned shape vectors and the principal 
components of the statistical face shape model to represent an approximate value of 

the test shape vector as x = x + £bj p x , where x is a mean shape vector defined 

j=i J 

according to the mean value of the feature points of the aligned shape vectors, p* is j- 

th principal component of the statistical face shape model, and bj is a corresponding 
projection coefficient; 

using a 2D scaled rigid transform algorithm to align the test shape vector 

represented by x = M(a,0)J^x + Sbj Pj J + t , where 0 , a and t are a rotating 

angle, a scaling factor and a shifting vector respectively; 

computing a normalized distance of the aligned test shape vectors by 

d(x)= ]sSJ ;and 

assigning one candidate feature point corresponding to one, having the 
smallest normalized distance, of the aligned test shape vectors as facial feature of the 
test face image. 

15. (Original) The method as claimed in claim 10, wherein the determining 
step of the second procedure further comprises the steps of: 
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using the average value of the aligned shape vectors and the principal 
components of the statistical face shape model to represent an approximate value of 

— k — 

the test shape vector as x = x + SKp; > where x is a mean shape vector defined 

j=i 

according to the mean value of the feature points of the aligned shape vectors, p* is j- 

th principal component of the statistical face shape model, and bj is a corresponding 
projection coefficient; 

using a 2D scaled rigid transform algorithm to align the test shape vector 

represented by x = M(ct,9)|^x + Zbj p- J + t , where 9 , a and t are a rotating 

angle, a scaling factor and a shifting vector respectively; 

computing an error value between the test shape vector and the average value 
of the aligned test shape vectors by 

B(x) = w, £|| I, (X) - t, - tb)p] \\ 2 +w 2 d(x) , where ±\\ I. (x) - t s - £b)v\ || 2 
i=i j=i ;=i j=i 

is a similarity of the test shape vector to corresponding aligned shape vector x a , and 

d(x) is the normalized distance of the aligned test shape vectors; and 

assigning one candidate feature point corresponding to one, having the 

smallest error value, of the test shape vectors as facial feature of the test face image. 

16. (Original) The method as claimed in claim 15, wherein the error value is 

computed by an equation e(x) = w, 



X) 



- w 2 d(x) . 
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